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Abstract: We investigate complex versions of the Korteweg-deVries equations and an Ito 
type nonlinear system with two coupled nonlinear fields. We systematically construct ra- 
tional, trigonometric/hyperbolic, elliptic and soliton solutions for these models and focus 
in particular on physically feasible systems, that is those with real energies. The reality 
of the energy is usually attributed to different realisations of an antilincar symmetry, as 
for instance 'PT-symmetry. It is shown that the symmetry can be spontaneously broken 
in two alternative ways either by specific choices of the domain or by manipulating the 
parameters in the solutions of the model, thus leading to complex energies. Surprisingly 
the reality of the energies can be regained in some cases by a further breaking of the 
symmetry on the level of the Hamiltonian. In many examples some of the fixed points 
in the complex solution for the field undergo a Hopf bifurcation in the PT-symmetry 
breaking process. By employing several different variants of the symmetries we propose 
many classes of new invariant extensions of these models and study their properties. The 
reduction of some of these models yields complex quantum mechanical models previously 
studied. 



1. Introduction 

One can adopt various points of view with regard to the usefulness of the study of complex 
classical and quantum mechanical systems. Being very orthodox one may just view the 
complex systems as providing a larger setting which allows a better insight from a broader 
framework when restricted to the real physical system. A very successful example for this 
viewpoint is the more than seventy year old proposal of the analytical S-matrix 
which is still pursued nowadays; especially in 1+1 dimensions. Genuinely non-Hermitian 
systems such as dissipative ones are also well studied, but they are usually regarded as 
open and are therefore not self-consistent ||. In contrast, a more recent perspective allows 
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to regard certain complex quantum mechanical Hamiltonians also as perfectly acceptable 
self-consistent descriptions of physical systems , in the sense that they possess real energy 
eigenvalue spectra and well defined unitary time-evolution (see || ||, f| for recent reviews). 
To consider pseudo/quasi Hermitian systems provides a very clear conceptual view in this 
respect as it makes use of a similarity transformation towards a Hermitian Hamiltonian 
system for which everything is well defined in a standard conventional sense, although such 
transformations do not exist for all types of such systems. A more radical view is to give a 
direct meaning to the non-Hermitian Hamiltonians without any reference to the Hermitian 
system. This latter point of view has to be taken in order to explain recent experiments in 
which non-Hermitian systems are studied on optical lattices, see e.g. ||, as the observed 
gain and loss can not be explained in a purely Hermitian setting. A further experimental 
realization has recently been proposed for graphene nanoribbons, where non-Hermitian 
Hamiltonians arise as effective Hamiltonians ||. 

Largely inspired by the study of the quantum systems, also classical complex systems 
have been investigated recently. Naturally one may view them too in various ways, either 
as directly meaningful or only sensible when transformed to a real system. Many classical 
models have already been investigated from these various perspectives, such as complex 



extensions of standard one particle real quantum mechanical potentials [1C, [TT], 12, 13], 
non-Hamiltonian dynamical systems [Q, chaotic systems fl5| and deformations of many- 
particle systems such as Calogero- Moser-Sutherland models g|, 0' HE HI- These 

inves- 
tigations led to new insights into the quantum theories based on the features found in 
these classical models, such as tunneling behaviour [2^], band structures J2l| and even 
a complex generalizations of Bohr's correspondence principle has been formulated p^ j. 
Extensions of field theories of nonlinear wave type, such as Korteweg-deVries (KdV) 



equations 23, 24, 25] and closely related models [26, 27, 28, 29] have also been investigated. 
These type of models will be the main subject of our investigations in this manuscript. 

Pseudo or quasi-Hermiticity are often equivalent to a simultaneous parity and time 
reversal, so-called 'P7~-symmetry. Remarkably this property, or more general antilinear 
symmetry [p0[|, of quantum mechanical systems is already visible on the classical level. It 
is known for more than fifty years that when this symmetry is unbroken the reality of the 
spectrum of the theory is guaranteed. Unbroken refers here to the property that 

[VT,H}=0 and PT$ = $, (1.1) 

i.e. the Hamiltonian H is PT-symmetric and its eigenfunctions $ are also eigenstates 
of the VT or any other antilinear operator. In case the latter property does not hold, 
one speaks of spontaneously broken T'T-symmetry and the eigenenergies become complex 
conjugate pairs. The 'PT-symmetry can be broken spontaneously in two different ways, 
either by tuning some of the parameters of the models appropriately or by manipulating 
the domain on which the model is defined, see e.g. |H]] for the latter possibility. One speaks 
of broken PT-symmetry when neither of the relations in ( |l . l| ) hold. In this case one usually 
expects complex energies, but we will demonstrate in this manuscript that in some cases 
real energy spectra can be produced by breaking the spontaneously broken theory further 
in a controlled way. 
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Here we will trace these properties on the classical level and use the reality of the 
classical energy as a criterium to select physically meaningful complex models including 
their boundary and initial conditions. The role of $ in ( |1.1| ) is in this case played by 
the solution of the classical field equation of motion, say u(x,t). For a given Hamiltonian 
density H depending on u(x, t) the energy on the interval x G [—a, a] is computed by 



E= r U \u(x)\ dx = 1% [u{x)\ —. (1.2) 

J-a Jt U x 



Here and throughout the manuscript we use the standard convention u x = alu/alx. When 
the system possesses periodic solutions, that is when u(a) = u(—a) along a path T in the 
complex u-plane, we can employ the alternative contour integral version in ( |1.2| ). A simple 
argument [^4| shows that the energy in the interval x £ [—a, a] is guaranteed to be real 
if the symmetry property %* [u{x)\ = % [u(—x)] holds for the Hamiltonian density. We 
will present here some unexpected examples of real energy solutions for which neither of 



the relations in (LI) hold, i.e. the T^T-symmetry is broken for the Hamiltonian and the 
solutions. 

Here we will investigate two different types of complexified wave equations. First of all 
we study in section 2 the most immediate way to complexify wave equations by introduc- 
ing complex boundary values and initial conditions for the KdV-system. We investigate 
systematically the traveling wave and soliton solutions. We study complex trajectories 
in the complex plane of the KdV-field especially with regard to their properties under 
'PT-symmetry breaking. In section 3 we employ the different types of 'PT-symmetry as 
a construction principle to propose new extended versions of the KdV-system. We inves- 
tigate the new models in a similar fashion as their undeformed counterparts. Particular 
attention is paid to the question of what kind of conditions will lead to physical models, 
in the sense that they possess real energies. In section 4 we study the effect of complex 
boundary conditions and initial values on a nonlinear system with two fundamental fields 
coupled to each other, which is referred to as Ito system for some specific parameter choice. 
We also investigate the P7~-symmetry properties of these models. Due to the presence of 
an additional field when compared with the KdV-system we can identify four different 
versions of P7~-symmetry being realised in these systems, which we exploit in section 5 to 
construct new models. We study them from similar points of view as the previous ones. 
We draw our conclusions in section 6. 



2. Complex Korteweg-deVries equation 

Complex extensions of the KdV equation have been investigated already some time ago 



for instance in 33, 34, 35 1. However, PT-symmetry has only been utilized recently in 



1 23, 24, 25 1 in order to understand some of their properties and in particular to construct 
new models. The standard KdV system is known to be a Hamiltonian system with density 



%KdV 



^3,72 



13,1 € 



(2.1) 
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leading to the Korteweg-deVries equation in the form 



u t + f3uu x + ju xxx = 0. 



(2.2) 



Usually the constants (5 and 7 are chosen to be real, but here we allow them to take complex 
values, thus including the possibility for u(x, t) to be complex. "PT-symmetry may then 
be realised in two alternative ways as 



both possibilities guaranteeing that VT± : %KdV ^ %KdV holds. The underlying models 
respect one of the two symmetries and are therefore different as they correspond to two 
distinct choices of the coupling constant which may, however, be related by a simple rotation 
in u. Nonetheless, with regard to possible deformations to be discussed below this second 
symmetry allows to construct different types of new models. 

A crucial feature of the model to be acceptable as physically consistent is the reality of 
the energy ( |1.2[ ) . In case of the Hamiltonian density Ukav the reality of the energy (|1.2j ) 
is guaranteed for the two possibilities u*(x) = u(—x) when ^,7 6 1 or u*(x) = —u(—x) 
when i/3,7 G R, resulting from VT+ or "P7 - -, respectively. 

2.1 PT-symmetric, spontaneously broken and broken solutions 

Let us now see which complex boundary conditions and initial values are physically per- 
missible. In order to establish this, we first briefly recall how the traveling wave solutions 
of the KdV-equations of motion may be constructed systematically. Integrating fl2,2|) twice 
leads to the equation 



with integration constants Ki,K2 G C and P(u) denoting a third order polynomial in u 
multiplied by an overall constant A. A further integration yields 



where we made the usual assumption that the field u(x,t) acquires the form of a traveling 
wave u(x, t) = u(Q with C = x — ct and c denoting the wave speed. In complete generality 
this is an elliptic integral, but it is instructive to generate simpler solutions by systemati- 
cally making some specific assumptions on the factorization of the polynomial P(u). This 
solution method may then be extended to the deformed cases. 

Demanding specific boundary conditions will impose further restrictions or might be 
entirely incompatible with certain factorizations of P(u). For instance, in case we wish to 
implement vanishing asymptotic boundary conditions for u and its derivatives, the once 
integrated version of equation ( |2.2| ) and (^) implies 



VT+ : x h-> —x, 1 1 y —t, i 1 y —i, u \-t u for /3, 7 G K, 
VT- ■ x 1 — y —x, 1 1 — y —t, i 1 — y —i, u 1— > —u for if3, 7 G R 




(2.4) 




(2.5) 





lim u, ur = =^ k\ = K2 = 0. 

C->±oo ^ 



(2.7) 
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In the following we will study the solutions in the complex u-plane. For this purpose 
it is useful to separate u into its real and imaginary part u R and u 1 , respectively, and 
decouple ( |2.5| ) into two first order differential equations in these variables 



,R 



±Re 



y/XJ P(u R + iu 1 ) 



and 



± Im 



V\\ P(u R + iu 1 ) 



In this set up we may then apply many of the techniques which have been developed for two 
dimensional dynamical systems, see for instance ||36f| . Most immediate is the application 
of the linearisation theorem at some fixed point Uf, converting the nonlinear system into 



u 



J(u R ,u% 



u 



(2.9) 



with Jacobian matrix 



J(« fl ,« J )| 

V ' >\u=Uf 



dRe[^y/P{u)} 

± ft?? 

± ft?* 



dKe[V\y/p{u)} 
± ft? 

dIm[V\y/P(u)] 
± ft? 



(2.10) 



We denote the eigenvalues of J(u = uf) by ji, j%. Provided the system (2.9) is simple and 



Re jj the linearisation theorem applies, stating that the phase portraits of the systems 
J2.8|) and (^) are qualitatively the same in some neighbourhood of the fixed point Uf. 
The ten similarity classes for 2 x 2-matrices fully characterizing all possible behaviours for 
the fixed points of (2-8) are reported for reference in appendix A. 



2.1.1 Rational solutions 

Factorizing P(u) at first in the simplest way as P(u) = (u — A) 3 with one constant A leaves 
no freedom when solving ( |2.5| ), as all constants are fixed 



A 



37 



w K2 



6/3 2 



and 



(2.11) 



Clearly asymptotically vanishing boundary conditions Q2.7| ) are only possible for a static 
solution with c = 0. The evaluation of the expression ( |2.6| ) produces for this factorization 
the rational solution 



"(0 = ^ 



127 



(2.12) 



f> /?(c-c ) 2 ' 

with additional integration constant Co £ Taking C as purely imaginary maintains the 
PT-symmetry of the solution, whereas any real part may be compensated by a shift in £, 
which is kept to be real. Independently of the parameter choice, asymptotically all rational 
solutions of the type ( |2.12| ) end up at lim^-too it(C) = c//3, as we can also observe in figure 

Focussing at first on the PT-symmetric scenario, we notice that while keeping the 
speed of the wave c real there are two possible choices for the conjugation leading to 
physical solutions, namely u*(() = u(—Q when z(q,/3,7 6 1 or u*(Q = ~ U (~C) when 
i£o,i/3,7 G M. We depict some complex trajectories in the u-plane in figure |l]a for several 
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different initial conditions Co- We have taken the plus and minus sign in fl2.6|) for the upper 
and lower half plane, respectively. We observe that for a specific branch the point A in the 
u-plane appears to be either a stable or an unstable improper asymptotic fixed point. We 
adopt here the notion for the characterizations of fixed points from the linearisation (see 
appendix A for a classification), despite the fact that is not a meromorphic function and 
the system is not easily linearized. This implies that only the choice with different branches 
for the upper and lower half will give rise to closed orbits as depicted. For increasing £, 
they run either out of the fixed point in the upper half plane and into it in the lower half 
or vice versa. The crossing of the trajectories with the real line is easily computed to be 
at u(0) = c/j3 + 127//3(Im£ ) 2 . This makes it evident that only complex trajectories may 
be closed, whereas the real solution is the only trajectory drifting off to minus infinity. 

Since the rational solution ( |2.12| ) does not have any free parameter left, as specified 
m (EHl , there is no possibility to break the PT-symmetry spontaneously for this type of 
solution. We may, however, break the PT-symmetry completely directly on the level of 
the Hamiltonian by fully complexifying /3 or 7, an example of which is depicted in figure 
|l]b. As expected we observe that the symmetry u*(C) = u(— £) has been lost, but instead 
the trajectories are almost symmetric about the line passing the two points A and u(0), 
which are now both located away from the real axis. However, the nature of the fixed 
point, being either an unstable or stable improper node, has not changed. 




Figure 1: Complex rational solutions of the KdV equation for different values of purely complex 
initial conditions £ : (a) PT-symmetric solutions for c = 1, /3 = 2, 7 = 3 and A = 1/2; (b) Broken 
PT-symmetric solutions for c = 1, j3 = 2 + i2, 7 = 3 and A = (1 — i)/4. Different trajectories are 
characterized by different initial conditions £ . Some values for the imaginary part of Im£ < 1 
are indicated on the trajectories. The unresolved white region corresponds to values for Im£ > 1. 
Throughout this manuscript we label panels from the left to the right as (a), (b), (c), etc. 



Considering the expression for the energy ( |l.2| ) it is evident that it will be real when 
we have the symmetry u*(x) = u(—x). We also compute the expression explicitly for the 
solution ( 2.12| ) by substituting it into ( |1.2j ). Then the energy in the interval [—a, a] is 
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computed to 



ac 

W 2 



c + 



367 



C 2 



+ 



72 7 2 
15/3 2 



10c (a 3 + 3aCo) 4 §7 (a 5 + 10a 3 Co + 5aCo) 



(« 2 -Cg) 3 (« 2 -Cg) 5 

(2.13) 

Evidently E a is real even when iCcb/-^7 £ ^ or i£g,z/3, 7 G K and complex otherwise, that 
is for the 'PT-symmetric and broken T^T-symmetric case, respectively. 



2.1.2 Trigonometric/hyperbolic solutions 

As the next possibility for the factorization we specify P(u) 



Af 



B) involving 



now two constants A and B, thus leaving one of them at our disposal when solving ( |2.5| ) 



L 
3 7 ' 



Ki 



2c), K 2 



A 2 

— (3c 
6 V 



2j3A) and 5 



3 4-2A. 



(2.14) 



Having now some freedom in the choice of the constants one may ask which ones are the 
most natural to use for the symmetry breaking. As we will see the constants A and B 
have a direct physical meaning and it appears therefore natural to view them as the free 
parameters to tune for a concrete model with fixed coupling constants /3,7, rather than the 
integration constants K\ or k 2 emerging more indirectly without immediate interpretation. 
By ( |2.7| ) vanishing asymptotic boundary conditions require the choice ^4 = 0. 

In general, the solution to ( |2.6D produces in this case the trigonometric/hyperbolic 
solution 



u (C) = B + (A - B) tanh 2 



2 



Wa(c-Co) 



(2.15) 



Let us first discuss the PT-symmetric scenario for which all constants are taken to be real 
except for £ , which we still allow to be complex. When in that case either A < B, A > 
or A > B, A < we obtain a periodic solution with period T = 2ir/( x /\A — B\y/\X\) as 
depicted in figure ||a. The closed trajectories surround the point A, whereas the point 
B is situated on its outside, following from the fact that on the real axis we always have 
u(C) < B or u (C) > B, in the respective cases. This behaviour is also confirmed by the 



linearisation ( [2.9| ) at the fixed point A. Parameterizing A 
the Jacobian ( |2.10| ) is easily computed to 



B = r A Be idAB and A 



>t>\ 



J(u)| 



u=A 



__ ( ±^r AB r\ cos [1(6 ab + 0\)] Ty/rAB~r x sin [\ 
y ±^/r A Br\ sin [\(9ab + 0\)] ±^r AB r x cos [± 



. AH + Ox)} 

sJtabt\ cos [\{6 A b + Ox)] 



(2.16) 



with eigenvalues 



ji = ±^r A Br\ exp 



?ab + Ox) 



and j 2 = ±^r A Brx exp 



-AOab + Ox) 



(2.17) 



Clearly for the two cases here, that is A < B, A > or A > B, A < 0, the eigenvalues of 
the Jacobian are purely complex, i.e. ji/ 2 = ±i\/\A — ~B\<y\\\, indicating that the point 
A is a centre. The result is independent of the sign in ( |2,8D . The linearisation at the fixed 
point B does not exist, due to the occurrence of the square root. 
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For this periodic case we may also compute the energy Et for one period T from 
-tt/vTA - S[i/fA[ to n/y/\A - analytically to 



# T = <fn[u(0) — = I 




(2.18) 



Here the contour T is taken to be any complex trajectory resulting for i£ E R. Then 



(2.18) follows from Cauchy's residue theorem and the fact that B will always be outside 
the contour T. We find that the energy is real and takes the same value for all trajectories 
independently of the concrete choice for the initial condition £ - Notice that the energy for 
the real solution, i.e. Co ^ ^> can n °t be computed in such an easy way as in that case T 
does not form a closed contour. Thus demanding PT-symmetry leads inevitably to purely 
complex initial value conditions and constitutes a natural "e-prescription" to deform the 
real solution with the purpose to compute the energy for one period. 

For the remaining possibilities ^4>i?,A>0oryl<i3,A<0ail solutions tend 
asymptotically to A, which we depict for some examples in figure |2|b. The linearisation 
Q2.9p at the fixed point A yields in these cases the two degenerate real eigenvalues j\ = j% = 
\/\A — B\ v|A| or j\ = j'2 = — y\A — B\ y/\X\ for J depending on the plus or minus sign 
in (|2.8| ), respectively. The Jacobian ( 2.16|) is diagonalisable in these cases and therefore A 
is an unstable or stable star node. As for the rational solutions, this implies that only the 
choice with different branches for the upper and lower half will give rise to closed orbits 
as seen in figure |||b. Now the energies can not be computed in a simple manner as for the 
periodic case since the singularity at A is situated on the contour. 



E 




4 
Reu 



0.2 




3.0 
Reu 



Figure 2: Complex trigonometric/hyperbolic solutions of the KdV equation: (a) "PT-symmetric 
periodic solutions with c = 1, /3 = 3/10, 7 = 3, A = 4, B — 2 and T — 2vT57r; (b) PT-symmetric 
asymptotically constant solutions with c = 1, /3 = 3/10, 7 = —3, A = 4 and B = 2. 



Let us next embark on the case in which the PT-symmetry is spontaneously broken, 
which unlike as for the rational solutions, is possible for ( [2.15 ) since we have additional 



- 8 - 



PT in complex nonlinear wave equations 



parameters at our disposal. From ( 2.14| ) and ( |2.15 ) it is clear that when 2(3 k\ + c 2 < the 
constant A becomes complex and we no longer have the property u*(C) = U (~C) ensuring 
the reality of the spectrum. The nature of the fixed point is now changed again to an 
unstable or stable focus depending on whether Rejj > or Re i; < 0, respectively. Thus 
once again closed orbits are obtained with the choice of different branches for the upper 
and lower half plane. This means of course that the periodic solution ceases to be periodic. 

The energies Et for one period in the periodic solution corresponding to the two choices 
A,B,(o and A*,B*,(o occur in complex conjugate pairs. This is the typical scenario for 
spontaneously broken PT-symmetry, i.e. the Hamiltonian is still PT-symmetric but the 
solutions to the Schrodinger equation in the quantum case and in the classical case to 
the equation of motion, are not. Formula ( p. 18 ) may still be used in this case for the 
computation of the energy even though the singularities are moved away from the real axis 
as they still lie within the contour. Once again the result does not depend on the choice of 
Co- We depict some trajectories for this type of spontaneous symmetry breaking in figure 
i The conjugate solution is simply obtained by a reflection about the real axis, as shown 
explicitly for the asymptotically constant solution in panel (b). 



E 




E 0.0 




Figure 3: Complex trigonometric/hyperbolic solutions of the KdV equation: (a) Spontaneously 
broken PT-symmetry of the periodic solution with c = 1, j3 — 3/10, j ~ 3, A = A + i/2 and 
B = 2 — i for Im£ = 0.5 solid (black), Im£ = 0.3 dotted (green) Im£ = 0.1 dashed (blue); (b) 
Spontaneously broken PT-symmetry of the asymptotically constant solution with c = 1, /3 = 3/10, 
7 = -3 for A = 4-i/2, B = 2 + i ImC = -0.5 solid (black)and A = 4 + i/2, B = 2 — i, ImC = 0.5 
dashed (blue). 



Finally, there is of course also the possibility to break "PT-symmetry directly for the 
Hamiltonian itself. For instance when /3 ^ R and/or 7 ^ R we expect the energies to be 
complex. In that case the complex conjugate energy would be obtained by considering a 
new type of Hamiltonian with (3*, 7* and thus it does not arise from within one specific 
model. We depict an example trajectory for this scenario in figure 

We observe from figure ||]a that the periodic broken solution near the fixed point is 
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qualitatively the same as the one for the spontaneously broken case, namely a stable or 
unstable focus. This behaviour follows from the eigenvalues ( 2, 17] ), which indicate that 
there is no distinction at this point in whether the complexification of the ji result from a 
spontaneous or a complete breaking of the PT-symmetry. Note that the trajectories still 
close even though this is not shown in figure ||a, but this may be seen on a larger scaled 
plot. We find a similar behaviour for the fixed point of the broken asymptotic solution as 
depicted in figure [|b. 



E 




a o.os 



—0,10 hi 
20 1.5 




Figure 4: Complex trigonometric/hyperbolic solutions of the KdV equation: (a) Broken VT- 
symmetry of the periodic solution with A = 4, B = 2, c = 1, /3 = 3/10, 7 = 3 + i/2 and lm£ = 6; 
(b) Broken 'PT-symmetry of the asymptotically constant solution with A = 4, B = 2, c = 1, 
P = 3/10, 7 = 3 + i/2 and ImC = 1/2. 

So far we have simply broken the symmetry by choosing some more or less random 
complex value for 7. However, we can also carry out this process in a more controlled 
fashion producing periodic motion and even some non-Hermitian non- "PT-symmetry cases 
with real energies. First of all we observe from eigenvalues of the Jacobian ( 2.1 0| ) that the 
fixed point A becomes a centre when Bab + #a = which is also compatible with the 
solution ( 2.15 ) from which we notice that a periodic motion occurs when X(A — B) < 0. 
This allows for complex values of the parameters /?, 7, A and B. Combining this constraint 
with the last equation in ( 2.14| ) leads to 



sin < 



8111 ( 



■ exp (iO A ) ■ 



(2.19) 



h-Op- 0a) 

Thus, for a given model, that is for any fixed values of f3 and 7 we obtain a periodic motion 
around the point A, given by the expression in ( 2.19j ) for any value of 9a- Indeed this is 
the case as we observe for instance in figure ^a. When using the constraint Oab + S\ = tt 
for the periodic motion in the expression for the energy ( 2.1§| ) we find 

7T PA 3 



Ex 



B\ 



(2.20) 
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Hence, demanding that the energy is to be real leads to the further constraint 39a + $t 



0,7r. Implementing this in (2.19) yields 



for A 



sin 9~ 



\P\ sin (Oy- 20^/3) 



cxp 



(2.21) 



This means for a given model, that is for any fixed values of /3 and 7 we can find a 
point A, given by the expression in ( |2,21| ), around which the trajectory is periodic and 
the corresponding energy is real. This holds irrespective of whether f3 and 7 are real or 
complex, or in other words whether the PT-symmetry is intact or completely broken. In 
figure ||b we depict an example solution of ( [2.21 ) corresponding to a periodic motion with 
completely broken PT-symmetry but real energy. An obvious question to ask at this point 
is whether by breaking this symmetry the system has acquired a new kind of anti-linear 
symmetry, which could serve to explain the reality of the energy. A systematic study of 
this more general issue will be presented elsewhere. 





Figure 5: Complex trigonometric/hyperbolic solutions of the KdV equation with completely bro- 
ken 'PT-symmetry: (a) Periodic solution with complex energy Et = —10.518 + il.666 for c = 1, 
P = l/2exp(i7r/4), 7 = l/3exp(«7r/3), A = V3 sin(ll7r/60) exp(-wr/10) w 3.025 - i0.983 and 
B = 6exp(-ivr/4) - 2>/3sin(ll7r/60) exp(-nr/10) w -1.806 - i2.277; (b) Periodic solution with 
real energy Et = — 47r for c = 1, (3 = 1/2 exp(i7r/4), 7 = l/3expi7r/3, A = 2V3 exp(— in/ 12) s» 
3.346 - i0.896 and B = —^6(1 +i) « -2.449 - i2.449. 



In general we observe from the various cases studied in this subsection that when 
we let the parameters vary one of the fixed points for the periodic solution undergoes a 
Hopf bifurcation, i.e. it changes its behaviour from a stable focus to a centre and then 
to an unstable focus. As we can see from ( 2.1 7| ) this behaviour is governed by the sign 
of ±cos [^(#ab + #a)] this bifurcation can be achieved either by a spontaneous symmetry 
breaking, that is by varying the free parameter A, or by a complete breaking of the sym- 
metry, by varying /3 or 7. Both cases pass through the PT-symmetric realised by 
the centre. 
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2.1.3 Elliptic solutions 

Next we specify P(u) = (u — A){u — B){u — C) with three constants A, B and C leaving 
two constants at our disposal when solving ( |2.5| ) 

A = -J^, ki = \ [(3(A 2 + AC + C 2 ) - 3c(A - C)\ , (2.22) 

AC 3c 
n 2 = — [3c-P(A + C)] and B = --{A + C). (2.23) 

6 p 



Vanishing asymptotic boundary conditions reduce this case to the previous one as by fl2.7| ) 
they require either A = B = 0, A = C = 0oiB = C = 0. The evaluation of fl2.6| ) yields 
in this case the elliptic solution 



u (0 = A + (B - A) ns ; 



IV^4VA(C-Co) 



A-C 



A-B 



(2.24) 



with ns(z|m) being one of the Jacobi elliptic functions. Therefore u (£) is a double periodic 
function 

u(C)=«(C + wi+w 2 ) (2-25) 

with periods 

8 „fA-C\ , 16 T /C-B\ 
wi = -: and w 2 = » ; -j • 2.26 

/b^Av^ \a-bJ Vb^aVx \a-bJ 

Here K(m) denotes a complete elliptic integral of the first kind. As expected we recover 
the previous case when two of the constants A, B, C coincide as in these cases the periods 
become dependent on each other, i.e. lira a^b ^1/^2 = l/2sign(C — A), liuiA^c ^1/^2 = 
—zoo and Huib-^c ^1/^2 = 0. 

Studying at first the "PT-symmetric solutions, we find periodic solutions encircling 
some of the fixed points. For instance in figure || we depict an example in which the fixed 
points A and B are encircled, whereas C is situated outside of the trajectories. Different 
types of scenarios are also expected, changing for instance the sign in 7 will lead to a 
periodic orbit surrounding the points B and C. 

Next we break the "PT-symmetry spontaneously by complexfying the free parameters 
of the solutions. We find that the periodic trajectories open up and cease to be periodic. 
In figure |7|a we trace part of the trajectory to illustrate the behaviour. Starting at the 
point u(— 64) ~ 7.19 + z0.74 the trajectory passes down between the points A and B and 
moves then up again surrounding once the points C and A in a clockwise sense. Thereafter 
it encircles C once more but instead of moving around A it passes inbetween C and A 
surrounding A in an anti-clockwise sense. It keeps progressing in an anti-clockwise manner 
encircling C before passing from below between the points A and B reaching the point 
u(18) ~ 6.36 + i3.45. It appears that this type of movement is repeated indefinitely. In 
figure 0b we depict a wider range for £ indicating that the region of the phase space depicted 
will eventually be filled densely by the trajectory, hence suggesting a chaotic behaviour. 
However, we do not observe a sensitivity towards the initial condition, which would be 
typical for a proper chaotic system. 
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Figure 6: 'PT-symmetric complex elliptic solutions of the KdV equation with A = 1, B = 3, 
C = 6, c = 1, p = 3/10, 7 = —3 for different values of Im£ . 




Re ii Re u 



Figure 7: Spontaneously broken PT-symmetric complex elliptic solutions of the KdV equation 
for ImCo = 6 with A = 4, B = 5 - i/2, C = 1 + i/2, c=l,0 = 3/10 and 7 = 3: (a) -64 < £ < 18 
solid (red) and 18 < £ < 200 dashed (black); (b) -200 < £ < 1400. 

Finally we may also break the 'P7~-symmetry completely by complexfying the param- 
eters of the model (3 or/and 7. Examples for such a scenario are depicted in figure ||. The 
behaviour is similar to the one of the spontaneously broken case, i.e. the periodic motion 
has turned into open trajectories with a noncompact limit set. Increasing the range for £ 
will fill the depicted part of the phase space similarly as in the spontaneously broken VT- 
symmetry case, thus suggesting a chaotic behaviour, albeit once again we do not observe 
the typical sensitivity towards the initial conditions. 

Figures ^ and || are very reminiscent of the plots which may be found in section 5 of 
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Figure 8: Broken 'PT-symmetric complex elliptic solutions of the KdV equation for Im £ = 6: (a) 
A = 1, B = 3, C = 6, c = 1, $ = 3/10 and 7 = 3 + 2i for -200 < C < 200;(b) A = 1, B = 2 + 3i, 
C = 6, c = 1, p = 3/10 - i/10 and 7 = 3 for -200 < C < 200. 



reference [12]. This is not surprising as formally the differential equations solved in there 
for the quantum mechanical setting for the potential V ~ x 3 are special cases of our more 
general treatment when making the identification u — )■ x and £ — > t for our traveling wave 
equation. With the further identifications 

n\ =0, k>2 = jE, j3 = Gcg and 7 = — c (2.27) 

equation (|2.5| ) converts precisely into the quantum mechanical Hamiltonian 

H = E = \p 2 + \x 2 -gx\ (2.28) 

treated in . The identification ( [2.271 ) explains why no analogue to our rational solution 
was found in [jT^I , since k\ = implies c = and therefore the vanishing of K2, i.e. the 
confining potential, and A. However, as can be seen from fl2.14j) there should be an analogue 
to our trigonometric solution with energy E = -4c 3 /7/3 2 when setting A = 2c//3. This 
energy depends on the coupling constant and is of course not freely choosable as in the 
elliptic case presented in this section where K2 is a free parameter. More potentials of the 
type treated in [jl2[ can be obtained systematically from the deformations discussed below. 

The new elliptic solutions recently found for generalized shallow wave equations [37] 
are expected to exhibit a similar behaviour as the solution discussed in this subsection. 

2.1.4 Soliton solutions 

There exist various techniques to construct soliton solutions. Here we briefly recall and 
apply Hirota's direct method |58] for that purpose. The starting point of this approach is 
Hirota's bilinear form, which reads in general 

p (D)t-t = 0, (2.29) 
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with p(D) being a polynomial in the Hirota derivatives, i.e. D x f ■ g = f x ■ g — f ■ g x where 
the dot "■" indicates the noncommutative nature of the product / • g. The soliton solutions 
to ( |2.29| ) may then be constructed perturbatively order by order in e from the expansion 

oo 

T = J2e n r n . (2.30) 

n=0 

Defining the function r(x,t) via the relation 

127 

u (x,t) = — (hiT) xx , (2.31) 



the KdV-equation (2^2) can be brought into Hirota's bilinear form 

( ^( 1 Dt + D x D t )T-T = 0. (2.32) 

We have set here the integration constant to zero. Taking the zeroth order solution to be 
To = 1 the first order tau-function results to T\(x, t) = exp(pix — ^p\t + a ) and we obtain 
the well-known one-soliton solution up to this order 

u (x, t) = 5 r1 - 3lPl = t. (2.33) 

/3cosh l2(Pix - JPp + 4>i)\ 

A XT-symmetric variant of this solution is depicted in figure ^a. As expected the origin 
in the u-plane is an asymptotic fixed point, reached for x — > ±00 or t — > ±00. 

Breaking the XT-symmetry in various ways yields more unexpected results. Regarding 
Pi) 4>ii Co as parameters of the solution and f3, 7 as model defining constants, their com- 
plexification constitutes a spontaneous or complete XT-symmetry breaking, respectively. 
The corresponding solutions in the it-plane appear as distorted versions of the symmetric 
case, as can be seen in figure [^b. An interesting feature of the broken complex solution 
is that, unlike its real counterpart, the soliton does not maintain its overall shape while 
traveling. Instead, it becomes a breather and regains its shape after a certain distance A x 
and certain amount of time At governed by the equation 

u(x + A x ,t) = u(x,t + A t ). (2.34) 

Separating 7, p\ into their real and imaginary part as 7 = -y r + i r ^ i and p\ = p r + ipi, we 
can solve ( 2.34| ) for the one-soliton solution ( p. 33 ) to 



At ~~ T~i J\ TT^ — m - — T\ — ' (2.35) 

[Pi -Prjli- 2 PiPr {Pi + Pr-)lr 

A o Pi far ~ Pi) 7i + 27T Pr fat ~ Pr) Jr (0 

[Pi ~ Pr) 7» - 2 PiPr {Pi + Pi) 7r 



The speed of the soliton is therefore given by 



£ = (3p? - rf ) 7r -«MzMh£. (2 . 37 ) 

t-±t Pr 
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Figure 9: Complex one-soliton solution of the KdV-equation at fixed time t = —2: (a) VT- 
symmetric solution with /3 = 6, 7 = 1, p% = 1.2 for <fi = i0.3 solid (blue), <fr — i0.8 dashed (red) and 
(\> = il.l dotted (black); (b) Broken 'PT-symmetric solution with /3 = 6, 7 = 1 + i0.4, pi = 1.2 for 
4> = i0.3 solid (blue), <f> — z0.8 dashed (red) and 4> — il.l dotted (black). 



Clearly for the T^T-symmetric solution the shape will be the same in the u-plane even 
at different times, whereas the complex solution will change its shape as we observe in 
figure ^ In figure 10 we observe that the one-soliton changes its shape while progressing 
in time but regains it after At = —tt/2 for the chosen values. 
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Figure 10: Complex spontaneously broken one-soliton solution of the KdV-equation with (3 = 6, 
7 = 1 + i/2, pi = 2, <fi = z0.8 and A t = — tt/2 for different times t — —ir/2 solid (blue), t = — 1 
dashed (red), t = dasheddot (orange), t — 0.7 dotted (green), and t = tt/2 dasheddotdot (black) 
(a) real part; (b) imaginary part. 



Proceeding to the next order in e we compute the two soliton solution to 

247ELo^(-l)W- fe 



u (x, t) 



P(pi +P2) 4 2 cosh [\ (r] l - rj 2 j) + e 2 



(P1+P2) 4 



+ 1 



(2.38) 



where we abbreviated r\ i = pix — jpft + (p i for i = 1, 2 with 



Co = 1 + cosh n 2 , ci=4sinhr/ 2 , C2 = cosh rj 1 + 6 cosh rj 2 — 1, C3 = 4 (sinhry 1 + sinh 77 2 ) 

(2.39) 

and used the symmetry relations Ci(7] 1 ,r] 2 ) = c^-.i(rj 2 .,r] l ). We depict the time evolution of 
this solution in figures 11 and |j~2|. 
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Figure 11: "PT-symmetric two-soliton solutions of the KdV equation for /3 = 6, 7 = 1, pi = 1.2, 
p 2 = 2.2, (f> 1 = iO.l and 4> 2 = i0.2. (a) t = -2 solid (blue), t = -0.2 dashed (red), t = 0.2 dotted 
(black); (b) t = 0.3 dotted (black), t = 0.8 dashed (red), t = 2.0 solid (blue). 



We observe in figure 11 the typical scenario for soliton scattering, albeit in the complex 
plane. For large negative time the two-solitons are separated, indicated here by two indi- 
vidual one-soliton solutions resembling the type depicted in figure || touching each other 
only in the asymptotic point at the origin. In the scattering regime the two solutions merge 
in a non T^T-symmetric manner until they separate again for large positive time, acquiring 
once again their individual one-soliton shape. 




Figure 12: Broken 'PT-symmetric two-soliton solutions of the KdV equation for f3 = 6, 7 = 
1 + in/8, p 1 = 2(2/3) 1 / 3 , p 2 = 2,<p l= iO.l and 2 = i0.2. (a) t = -4 solid (blue), t = -3.5 dashed 
(red), t = -2. dotted (black); (b) t = 0.7 solid (blue), t = 2 dashed (red), t = 8 dotted (black). 



For the values in figure ^ the periods for the one-soliton breather to regain its shape 
with pi and P2 are computed by ( 2.35| ) to be = — 3 and = —2, respectively. Indeed 
we observe in figure Il2|a that one of the solitons has recovered its shape comparing the 



solutions at t = —4 and t = —2. In figure 12b we observed that both solitons have almost 
regained their original shape passing from t = 2 to t = 8, that is after 2A| = 3A^. The 
slight discrepancy visible in the figure is due to the well known time delay occurring in a 
multi-soliton scattering, see e.g. |^|. 

The energy for the one-soliton solution on the interval (—00,00) is easily computed 
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Eis = -^P, (2.40) 



with formula (|l.^ 

367 3 pf 

exhibiting the typical behaviour, that is being real for the PT-symmetric case and complex 
for the (spontaneously) broken scenario. The integral ( |1.2| ) with u(x, t) given by two-soliton 
solution Q2.38 ) is not easily computed analytically. For large values of time we compute 



numerically real energies for PT-symmetric solutions, for instance for the values used in 



figure 11 we obtain Ei s ~ —10.8049 = E\ s {j>\) + E\ s {j>2). On the other hand we compute 



complex energies for broken PT-symmetric solutions, such as for the values used in figure 



12 for which we evaluate E2 S ~ —7.8876 — z9. 4327 = Ei s (pi) + Ei s {p2). This means that in 



all scenarios the energy of the two-soliton equals the sum of the two separate one-solitons. 
3. PT-symmetric deformations of the KdV equation 

Employing the standard arguments used in the study of non-Hermitian Hamiltonian sys- 
tems with real eigenvalues, we may now PT-symmetrically deform the Hamiltonian and 
maintain the possibility to have well defined physical systems, e.g. we still obtain real 
energies resulting from the new models despite their non-Hermitian nature. The general 
principle is simply to deform P7~-anti-symmetric quantities, that means if we have the 
property VT : 4>(x,t) h-> —(p(x,t) for some field 4>(x,t), we define a deformation map 5 E : 
<j)(x,t) i — y —i[i(j)(x,t)] £ . The undeformed case is recovered in the limit e = 1. The new 
quantity will remain anti-PT-symmetric with the crucial difference that the overall minus 
sign is generated from the antilinear nature of the PT-operator, i.e. i i— > —i, rather than 
from (j)(x,t) i — y —cf)(x,t). This means for the Hamiltonian resulting from the density ( |2.1| ) 
we can make use of either of the two possibilities 

'. Ux ' y u X: e := —i(iu x ) £ or 5~ : u h > u £ := —i(iu) £ , (3-1) 

depending on whether we choose u(x, t) to be PT-symmetric or PT-anti-symmetric, re- 
spectively. Accordingly we define the deformed models with some suitable normalisation 
by the following Hamiltonian densities 

H + = -^- T L-(iu x r+\ and = {l+ P {2 + e) (iuy +2 + 1 -ul (3.2) 

with corresponding equations of motion 

u t + (3uu x + ^u XX x,e = and u t + i(3u e u x + ju XX x = 0, (3.3) 

respectively. The Hamiltonian %f was proposed in [24] whereas % E corresponds to a 
complex version of the generalized KdV equation. For the higher deformed derivatives we 
use here the notation u xx ,e ■= d x u x , e , u XX x,e ■= d x u XjE ,- ■ ■ ,u nx ,e ■= d x ~ l u X)E , which means 
we only deform the first derivative and keep acting on it with d x to define the higher order 
derivatives. On the level of the equation of motion we could also deform the dispersion term 



proportional to /3, as investigated in |23J for the KdV-equation. However, such deformations 
do not lead to Hamiltonian systems and the question of how PT-symmetry may be utilized 
in this scenario remains an open issue. 

Let us now present some solutions to these equations. 
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3.1 The %+-models 



Proceeding similarly as in the previous section, we integrate the first equation in ( |3.3| ) twice 
and obtain the deformed version of equation (|2.5| ) 



cxp 



ITT 



2(8+1) 



(1 - e + 4n) 



[A eJ P(«)]- 



(3.4) 



where we abbreviated A e = —(3(1 + e)/(6je) and denote different branches by n. The 
polynomial -P(u) is identical to the one introduced in ( |2.5[ ), which means we can employ 
the same factorization as in the previous section. Integrating once more yields 



c w _ Co = exp 



ITT 



2(8+1) 



1 + An) 



du- 



1 



[KP(u)\ 



(3.5) 



We may now proceed as before and specify further the form of P(u). As for the case 8 = 1, 
in some specific cases we succeed to compute the remaining integral and subsequently solve 
the resulting equation for it, thus obtaining u(C). However, even when this is not possible 
analytically we can still investigate (^T^) numerically for all cases by viewing £ as a function 
of the complex variable u and plotting the contours of Im[C(u)] = Co m the tt-plane, while 
taking special care about the different Riemann sheets labeled by n. One should distinguish 
here these type of Riemann sheets from those arising due to the technique employed in our 
solution procedure. Considering £ as a function of u, as we do in some intermediate steps, 



will introduce new branch cuts which are sometimes seen in our figures, e.g. figure 27b. 



However, we do not attribute any meaning to them in the it-plane. Genuine branch cuts 
can always be distinguished from the "technical" ones by the fact that they have to be 
connected to the fixed points which are branch points. 

3.1.1 Rational solutions 

We start with the same assumption as in the previous section, namely P(u) = (u — A) 3 , 
which will impose the same constraints ( 2.11| ) for the factorization. The case e = 2 is 
special since we may take the root 1/3 in this case. Integrating (|3.5| ) and solving the result 
for u we find 



u 



(n) 



(0 = 4 + exp 



(C - Co) 



2 2/3 V7 

For the remaining integer values e E Z\{2} we can compute also a particular solution 

' iir 5 + e \ ( 67s 



(3.6) 



«(C) 



+ exp 



2 2 



e-2 



£ + 1 



(l + e)— (e-2)— (C-Co)~ 



(3.7) 



Our 



However, in general we have to take care about the different branches present in ( |3.5| ). 
numerical findings for some specific cases are depicted in the figures below. 

We supplement our numerical analysis with the prediction of some analytical features. 
Of special interest are the lines approaching or leaving the fixed points radially. When the 
symmetry is unbroken they can be defined equivalently as the lines for which Co = 0- We 
compute them by noting first that on one hand any point u on the line radially crossing the 
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point u = A is characterized by a constant value for arg(u — A) =: 9q. On the other hand 
the change of u with respect to £ has to point into the same direction, i.e. arg(±-u^) = 9q. 
Hence the lines of real initial conditions are determined by solving 



arg(n^ n) — A) + 2nm. 



Employing (3.4) and the factorized form of P{u) we can solve (3.8), obtaining 



□ (n,m,±) 
; 



arg(n 

7T 



2(2 -e) 



(») _ a) 
[e - 1 + 4(m ■ 



(3.8) 



(3.9) 



n + me) - (1 + e)(l ± 1)] + 



arg A £ 



for the angle in which the trajectories with £ = enter the point A. 



In figure 13 we show the first two members of the sequel parameterizing the deformation 
parameter as e = —nj (n+1) with n = 1, 2, . . . We observe that the pictures resemble flowers 
with 4 + 6n petals. We may predict the number of petals analytically as they are equal 
to the number of Co = solutions computable by ( |3.9|) . For instance for the case n = 3 
the formula Q3.9| ) predicts the 16 values #o = (2^ — 1)tt/16 for the parameter choice c = 1, 
P = 2 and 7 = 3 with t = 1, . . . , 16. We recognize these values in figure |l3|p. Note that 
the wedge regions separate different solutions from each other as the point A is always an 
asymptotic fixed point. 





Figure 13: Complex 'PT-symmetric rational solutions of the deformed KdV equation with A 
1/2, c = 1, f3 = 2 and 7 = 3 for (a) ft+ and (b) H+ 2/3 . 



The trajectories appear to be qualitatively quite different for negative integer values 
of e. For instance in figure 14a and [X4] b we depict some trajectories for the models with 
e = —2 and e = —3, respectively. It appears that the fixed point A is more like a saddle 
point, i.e. we exhibit trajectories more of a hyperbolic nature running away to infinity 
rather than converging asymptotically to A as in figure |l^. The C = solutions are 
predicted once again correctly by ( [3.9D to be at 8q = (2£ — l)vr/8 and 6q = (2£ — l)ir/5. 
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Figure 14: Complex 'PT-symmetric rational solutions of the deformed KdV equation with A = 1/2, 
c = 1, (3 = 2 and 7 = 3 for (a) H± 2 and (b) Ht 3 . 



For positive rational values of e a complete trajectory extends over several different 
Riemann sheets. Figure 15a and 15b show the solutions and C , respectively. The 
angles for the £ = solutions are predicted correctly by Q3.9| ) to be #0 = 4£ir/5 for 
£ = 1, . . . , 5. A closed trajectory is obtained when passing the branch cut at —00 to 1/2 
from the upper half plane in panel (a) to the lower half plane in panel (b). In figure |l7| a a 
single trajectory is depicted for Im£o = 1. 





Figure 15: PT-symmetric rational solutions for Hy 3 : Different Riemann sheets for A = 1/2, 
c = 1, & = 2 and 7 = 3 (a) C (1) and (b) C (2) - 



For the broken "PT-symmetry we present the solution in figure 16- The branch cut 
at —00 — i/4 to (1 — i)/4 is passed from above in panel (a) to below in panel (b). This 
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is illustrated in figure 17b, where we depict just one single trajectory for Im£ = 1. The 



trajectories for the PT-symmetric and broken PT-symmetric case look qualitatively very 
similar, the major difference being that the fixed point has moved away from the real axis, 
thus leading to a loss of the symmetry. 





Figure 16: Broken "PT-symmetric rational solutions for HJ, 3 : Different Riemann sheets for 
A = (1 - i)/4, c = 1, /3 = 2 + 2i and 7 = 3 (a) C (1) and (b) C (2) - 




-2-10 1 2-2-10 1 2 



Re u Re u 

Figure 17: Single trajectory for the rational solution of 'H 1 / 3 with ImC = 1 : (a) PT-symmetric 
solutions for the values as specified in figure [IB], (b) Broken PT-symmetric solutions for the values 
as specified in figure [l6| . 

Potentially there are of course many more possible values for e to be considered. We 
conclude here just by presenting two more examples with e being an integer, as these 
are the most common deformations usually considered. For these values we obtain yet 
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Figure 18: Single trajectory for the rational solution for of with Im£ = 1: (a) PT-symmetric 
solutions for A = 1/2, c = 1, f3 = 2 and 7 = 3, (b) Broken 'PT-symmetric solutions for A = 
1/4(1 - i), c = 1, P = 2 + 2i and 7 = 3. 



another type of characteristics as more and more branches have to be taken into account 
for increasing e. In figure ^ and 19 we depict all branches for the trajectory with Im£ ( 
for e = 3 and e = 6, respectively. 



1 




Figure 19: Single trajectory for the rational solution for of T-Lq with ImC = 1: (a) PT-symmetric 
solutions for A = 1/2, c = 1, j3 = 2 and 7 = 3, (b) Broken PT-symmetric solutions for A = 
1/4(1 - i), c = 1, P = 2 + 2i and 7 = 3. 



We observe some intricate winding behaviour near the fixed point, which, however, is 
not asymptotic in these cases. As is apparent from the solution (|3.7| ) the trajectories diverge 
to infinity into various directions depending on the chosen Riemann sheet. Evidently the 
number of these asymptotes grows with increasing e. Breaking the PT-symmetry will only 
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distort the trajectories, giving rise to new directions of the asymptotes, but not changing 
their numbers. 

3.1.2 Trigonometric/hyperbolic solutions 

As in the non-deformed case we assume next the factorization P(u) = (u — A) 2 (u — B), 



which will impose once more the constraints ( 2.14 ). Also in this case we can compute ( |3.5| ) 
numerically. 

Again we supplement our numerical findings by some analytical results. We can predict 
once more the lines for which £ = explicitly following the arguments of the previous 
subsection, with the difference that we have to take the limit to either the point A or B. 
For the angle of lines near the point A we find 



-)(n,m,±) 
1 A 



7T 



1 - e + 4(n - m) - 4me] H — &rg(A — B)-\ — arg(±A e ) (3.10) 



2(e-l) L " ' ~ J ' e- 1 

whereas for the angle of lines near the point B we compute 



g(n,m,±) 



IT 



(e - 1 - 4n) + (1 + e)2irm + (e - 1) arg(B - A) - arg(±A £ 



(3.11) 



We depict some specific examples with different characteristics in the figures EQ to 24 





Figure 20: Complex PT-symmetric trigonometric/hyperbolic solutions of the deformed KdV equa- 
tion with A = 4, B = 2, c = 1, ft = 2 and 7 = 3 for (a) U± 1/2 and (b) U + _ 



-2/3" 



In figure 20 we recognize a similar characteristic behaviour as for the rational solution 



of the same model depicted in figure 13. Roughly the portrait 20 corresponds to 13 with 



the difference that the single flower centre in the rational case situated as A has been 
pulled apart to the two points A and B in the trigonometric one. The overall effect is 
that some of trajectories separated in the rational case in different wedges are joint in the 
trigonometric/hyperbolic case. To illustrate this we have extracted one single trajectory 
in figure Ell We observe a similar behaviour for the entire sequel parameterized by e = 
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Figure 21: Single trajectory for the complex ■pT-symmetric trigonometric/hyperbolic solutions 



of the deformed KdV equation for 'H_ 1 / 2 with the same values as specified in figure 20 for (a) 
Im£ = 1 and (b) Im( = — 1. 



A 




A 




Figure 22: Broken 'PT-symmetric trigonometric/hyperbolic solutions of the deformed KdV equa- 
tion 'H+ :| y 2 : (a) Spontaneously broken "PT-symmetry with A = A + i, B = 2 — 2z, c = 1, /3 = 3/10 
and 7 = 3; (b) broken PT-symmetry with A = 4, B = 2, c = 1, /? = 3/10 and 7 = 3 + i. 



—n/(n + 1) with n = 1,2,... Notice also that unlike as in the case e = 1 there is no 
distinction between a periodic and an asymptotically constant solution. 

Having enough free parameters available we can now break the symmetry for this 



solution also spontaneously or completely as illustrated in figure |22|. In both cases we 
observe that the amount of wedges remains unchanged, the symmetry about the real axis 
is lost and the winding around the fixed points becomes more intricate. In addition, 
trajectories from certain wedge regions connect in different ways, e.g. in figure 2C we 
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observe a connected trajectory in the two light shaded wedge region to the left of the 
vertical as can also be seen in figure Such type of trajectories do not exist in the 



broken case as we can observe in figure 22. We note that instead trajectories from the light 
shaded region in the upper half plane connect to a dark shaded region in the lower half 
plane to the left of the corresponding light shaded region in the unbroken case. 




Figure 23: Complex "PT-symmetric trigonometric/hyperbolic solutions of the deformed KdV equa- 
tion with A = 4, c = 1, P = 2 and 7 = 3 for (a) H± 2 and (b) Ht 3 - 



i 




Figure 24: Broken 'PT-symmetric trigonometric/hyperbolic solutions of the deformed KdV equa- 
tion for H^ 2 '- ( a ) Spontaneously broken PT-symmetry with A = 4 + i/2, B = 2 — i,c = 1, j3 = 3/10 
and 7 = 3; (b) broken PT-symmetry with A = 4, B = 2, c = 1, /3 = 3/10 and 7 = 3 + 4i. 



Comparing next the solutions for e G M we observe in figure 23 a a similar transfor- 
mation from the rational to the trigonometric/hyperbolic case as for the previous example. 
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Again we find that the rational solution resembles this solution for the same values of the 
deformation parameters, with the difference that the characteristic behaviour around the 
point A in the rational case has been distorted to the points A and B. 



The broken PT-symmetric case for this model is presented in figure 24. In panel (a) 
we depict the spontaneously broken PT-symmetric scenario for %^ 2 ' observing that the 
symmetry is lost since the trajectories rotate together with the two fixed points. The PT- 



symmetry is completely broken in the presentation in 24 b for the same model. Having kept 



the values for the fixed points, we note that in this case the orbits rotate around them. 

We have investigated also other models observing similar patterns, which we will how- 
ever not present here. 

3.1.3 Elliptic solutions 




Figure 25: Elliptic solutions for H± 2 - ( a ) PT-symmetric with A = 1, B = 3, C = 6, /3 = 3/10, 
7 = —3 and c = 1; (b) spontaneously broken PT-symmetry with A — 1 + i, £> = 3 — i, C = 6, 
f3 = 3/10, 7 = -3 and c = 1; (c) broken PT-symmetry with A = 1, B = 3, C = 6, /3 = 3/10, 
7 = —3 + 4i and c = 1. 




K=u. fell 



Figure 26: Elliptic solutions for U\,^. (a) PT-symmetric with A = 1, B = 3, C = 6, /3 = 3/10, 
7 = —3 and c = 1; (b) spontaneously broken PT-symmetry with A = 1 + i, B = 3 — i, C = 6, 
£ = 3/10, 7 = -3 and c = l. 
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Proceeding just as in the undeformed case we specify next P(u) = (u—A)(u—B)(u—C), 
which imposes the constraints ( 2.22| ) and ( 2. 23] ), thus leaving two constants free. We 
present here only the results for the cases e = —2 and e = —1/2 in the figures 25 and 26, 
respectively. 

In all three possible scenarios we observe a similar qualitative behaviour as for the 
trigonometric/hyperbolic solutions of the previous subsection with the difference that we 
have three instead of two fixed points. 

3.2 The %7-models 

Let us now turn to the second type of deformation. In order to construct the solutions 
we proceed in a similar manner as in the previous case. Integrating now twice the second 
equation in equation ( |3.3| ) we obtain 

2 



1 



where 



c 2 

k 2 + K\U + -u 



A 



l + e)(2 + e) 
2pi e 



u 



2+e 



--: XQ(u), 



(3.12) 



(3.13) 



7 (l + e )(2 + e) 

A crucial difference between W£ and H~ is that unlike as the polynomial P(u), which was 
of fixed order 3, the order of Q(u) depends on e, meaning that we have more and more 
possibilities to factorize it for growing e. For instance, for a given integer value n G No the 
factorization of Q(u) as 

Q{u) = {u- Ai) e+2 - n TT. (u- A l+1 ), (3.14) 

admits solutions provided n — 2 < e < n + 1 and e G N. This allows of course for yet 
another infinity of possibilities. 

When ki = K2 = we can find a closed solution valid for all e by integrating ( 3.1 2| ) 
and solving it for u 



«(0 



c(e + !)(£ + 2) 



1/6 



i e j3 



cosh 



Vce(C-Co) 



+ 1 



(3.15) 



The generic scenario does not yield such a simple answer. 
3.2.1 H 2 

This case is especially interesting as it corresponds to a complex version of the modified 
KdV-equation. Specifying ( 3.14 ) for instance as Q(u) = (u — A) 3 (u — B) we can factorize 
the polynomial in (|3.12|) with the choice 



2c 3/2 



Kl 



K 2 



C 



A = -j and B = -^=. (3.16) 



?>J^ 4/9' 3 y/=fi 

Notice that this fixes all the boundary conditions for a given model. Solving ( |3.12j ) then 
yields a rational solution for the second equation in (|3 

-2 

«(0 



c 2< 2 - 97 



/? 3 7 + 2cC 



(3.17) 
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As is well know for the real case one may construct solutions for the KdV-equation from 
those of the modified KdV-equation by means of a Miura transformation. We expect this 
also to hold for their complex versions. Indeed, using the transformation of the form 



«KdV (C) 




u 



(3.18) 



yields the rational solution of the KdV-equations ( |2. 12| ) from ( ft.17 ) when we identify 
Co = i\/3j/2c therein. 

Assuming in ( |3.14| ) instead Q(u) = {u — A) 2 (u — B)(u — C) we can factorize the 
polynomial in fl3.12|) with the constraints 



A 



f3C 2 (g - 5Cf3) + 9c (g - 3C/3) 
81/3 ! 

-~(B + C) and B 



H2 



C{2§-Cf3) (C/3 + $Y 
324/3 2 



3/3 



(3.19) 
(3.20) 



where we abbreviated ?9 := \/2y / —(3 (f3C 2 + 9c). Notice that in this case one constant 
remains free. The integration of ( |3 .12 ) yields in this case a trigonometric solution, which, 
using the same Miura transformation ( 3.18j ), may also be converted into a solution of the 
KdV-equations. 

3.2.2 HI 



In this case the polynomial of the right hand side of ( |3.12 ) is of sixth order. We present 
here just one very symmetric solution by assuming a factorization of the form Q{u) = 
u 2 (u 2 — B 2 )(u 2 — C 2 ), which is possible with the simple choice 



Ki = K2 = 0, B = iC 



and 



4 _ 15c 

T' 



(3.21) 



Thus we have made contact with the solution ( |3.15| ). Parameterizing B = r^e* 6 ^ and 



A = r\e l ® x the eigenvalues of the Jacobian when linearized about u = are easily computed 



to 



n 



±i 



r A r B exp 



:(40B + e x ) 



and 



h = =FV^A r B ex P 



(3.22) 



This means for the PT-symmetric solution we always obtain two real degenerate eigenvalues 
and therefore a star node at u = 0. For the positive square root in ( 3.12j ) with B = 
±(15c//3) 1//4 or with B = ±i(15c//3) 1//4 the node is stable or unstable, respectively. The 
stability property is reversed for the negative square root. Taking the branch cuts at 
Imu = 0, Reu = 0, Imn = Reu and Imu = — Reu into account we obtain closed curves. 
The four zeros of Q(u), i.e. ±(15c//3) 1/ ' 4 , ±i(15c//3) 1//4 , are surrounded by the trajectories. 
All these features can be seen in figure p7|a. We may also tune the coupling constants in 
such a way that 48b + @\ = 0, which by ( |3,22 ) implies that the fixed point at the origin 
should become a centre. Indeed, for a specific choice we observe this in figure p7|b. Apart 
from the origin, the remaining four fixed points are now situated outside of the trajectories. 
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-4-5 3 4-4-5 3 4 



feu. fcu. 

Figure 27: 'PT-symmetric solutions for : (a) Star node at the origin for c = 1, j3 = 2, 7 = 1 
and B = (15/2) 1 / 4 ; (b) centre at the origin for c = 1, j3 = 1, 7 = -1 and S = (15/2) 1 / 4 . 




-4-3 3 4-4-3 3 4 



Figure 28: Broken PT-symmetric solution for H^: (a) Star node at the origin for c = 1, j3 = 2 + i3, 
7 = 1 and £? = (15/2 + J3) 1 ' 4 ; (b) centre at the origin for c = 1, (5 = 2 + z3, 7 = — 1 and 
B = (30/13 -^S/IS) 1 / 4 . 



Since we have no free parameter left in our factorization this solution can not be broken 
spontaneously. A complete breaking is carried out by complexifying f3 or 7. When choosing 
Im j3 7^ we can still achieve in ( 3.22[ ) that 46b+0\ = ±7r, such that the star node nature of 
the fixed point is preserved despite the fact that the 'PT-symmetry is lost. For a particular 
choice this behaviour is depicted in figure 28 a. More surprising is the fact that we can also 
achieve that 48b + #a = in the broken case. This means the eigenvalues in ( 3.22[ ) are 
purely imaginary and the fixed point at the origin is a centre. We depict this possibility in 
figure 28b. This means we have closed trajectories even in the P7"-symmetrically broken 
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case. 



In contrast for Im7 7^ the eigenvalues ji will become complex and the fixed point 
at the origin turns into an unstable or stable focus. Taking the branch cut structure into 



account we observed this behaviour in figure 2£a for j3 G M, 7 E C and in figure p9p for 
0€C, 7GC. 



B 





Figure 29: Broken 'PT-symmetric solution for % 4 : (a) focus at the origin for c = 1, j3 = 2, 
7 = l + i3 and B = (15/2) 1 / 4 ; (b) focus at the origin for ImC = 0.8, c = 1, /3 = 2 + i3, 7 = -l + i2 
and S = (30/13 - iAb/13) 1 ^ 4 . 

Once again we observed that when identifying a coupling constant as a bifurcation 
parameter, the origin undergoes a Hopf bifurcation. 

Making once more the identification u — > x and £ — > t for our traveling wave equation 
together with the identification 

ki=0, K 2 = <yE, and p = 75(1 + e)(2 + e) (3.23) 

equation ( 3.12| ) converts into the classical deformation of the harmonic oscillator 

(3.24) 



1 c 

H = E = -p 2 - —x 2 + gx 2 (ix) e . 



2 7 

Setting furthermore c = 0, which in our setting corresponds to a static solution, we obtain 
precisely the potential treated in the seminal paper by Bender and Boettcher 0]. 



4. Complex coupled nonlinear wave equation of Ito type 

It is known for some time that the Korteweg-deVries field u(x, t) may be coupled to a 
nonlinear field v(x,t) without destroying the integrability of the newly constructed system. 
This coupling is carried out by the introduction of an additional dispersion term involving 
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the second field v(x,t) and the assumption of a specific form of the evolution equation for 
this field. The coupled system describing this type of scenario acquires the general form 

u t + avv x + (3uu x + <yu xxx = 0, a,/3,jeC, (4.1) 

v t + 8{uv) x + 4>v xxx = 0, 8, <j> G C. (4.2) 

Imposing the constraints (3 = 37, (j) = 1, 8 = 3 on the constants and leaving a, 7 free, 
the Hirota method can be applied to establish that the system ( |4.l| ) and (|4.2| ) possesses 
N-soliton solutions Eo|. When selecting in addition 7 = —1/2 as in El] or a = —2, 



/3 = —6, 7 = —1, (5 = —2 and = as in [42|, it was established in both cases that the 



system possesses infinitely many charges. For = the system ( [4.1| ) and (4.2) was shown 



to possess soliton solutions of cusp type [43]. Notice also that when complexifying the 



KdV-field in (2.2) as uxdv u + iv we obtain the system ( f4.1| ) and ( f4.2j ) for the special 



choice a = -/3 KdV , f3 = /3 KdV , 7 = 7 K dV> 5 = /^KdV and 4> = TKdV- In the following we will 



mainly discuss the case eft = 0, such that even the real case of (4.1) and (4.2) is distinct 
from the complexified version of the KdV equation. 

It is straightforward to verify that for the choice S = a the system of equations ( f4.1| ) 
and ( [4.2[ ) results from a variation of a Hamiltonian whose density is given by 

1 = ~2 UV ~ 6 2 x 2 ^ ' 

when using the variational derivative and time evolution in the standard way 

^ = 9- x E n=0 (-ir^^-;) ** (4-4) 



The Hamiltonian resulting from the density (|4.3| ) is manifestly PT-symmetric as it remains 
invariant under a simultaneous parity transformation and time reversal which may be 
realised in four alternative ways as 

VT++ : x 1 — y —x, 1 1 — y —t, i h-> — i, u^u,i)4i; for a, j3, 7, G R, (4-5) 

PT-i : x 1 — ^ —x, i 1 — y —t, i 1 — y —i, m^ii,!)^ — v for a, /?, 7, </> G R, (4-6) 

"PT 1- : 2; h-> —re, 1 1—)- — i, i I—)- — i, u 1— )• — it, ui^d for za, if3, 7, G R, (4.7) 

PT : x I—)- —x, i 1 y —t, i 1 y —i, u \-t —u, v \-t —v for ia, ij3, 7, <f> G R, (4.8) 

depending on whether we choose the fields u, t> to be PT-symmetric or anti-symmetric. All 
possibilities ensure that VTij : Hi h-> holds with z,j G {+,—}• First we will exploit 
these possibilities to explain the reality of the energies and in section 5 we use them to 
define new physically feasible models. 

4.1 "PT-symmetric, spontaneously broken and broken solutions 

As usual in this context we assume that the fields acquire the form of a traveling wave 



u(x,t) = u(C) and v(x,t) = v(() with £ = x — ct. In [44| it was shown for the case 



that the possibility to have only one field to be a traveling wave and not the other is 
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inconsistent. We extrapolate here without rigorous proof that this assumption can be made 
without loss of generality even in other cases. Let us briefly recall how these equations 
may be solved in a systematic way. To begin with we integrate (4.1) and ( |4.2| ) to 



-cu + ^v 2 + ^u 2 + 7it££ = Ki, (4.9) 
—cv + auv = K2, (4-10) 

with integration constants k±,k,2- In the following we always exclude the case k<i = as 
this implies the vanishing of the new field v = 0, which means a reduction to the deformed 
KdV equation, or a constant KdV- field u = a/c. Multiplying ( |4.9| ) by u^, using ( |4. 10 ) to 
replace v by u in (14. 91), we can integrate once more and obtain 



u 2 c = -Ik 3 + k 1 u + -u 2 -^u 3 + ^ . 4.11 

7 V 2 6 2 au — c J 

This equation is difficult to solve directly, but following we trade the u field for the v 
field with 

u = 1 and ur = ^vt (4-12) 

a av av z 



and obtain 



where 



ao = Bk^, a\ = 3ck 2 (B — a), 02 = 3k2(/3c 2 — 2ac 2 — 2a 2 K,\), 

03 = c 3 (/3 — 3a) — 6a 2 (cKi + 0/^3), 04 = — 3a 3 «2- 



(4.14) 



At first sight this looks even less encouraging than ( (4.1l|) . However, now the right hand 
side is a polynomial and in case we can factorize the sum X^t=o a k vk some convenient 
form we may be able to integrate (|4.13j ) similarly as in the previous section. Up to one 
integration the solution is therefore 

±VA(C-Co)= / dv -±=, (4.15) 

with A = a 2 /jK2 arid XR(v) corresponding to the right hand side of (|4.13|) . We note that 
not all conceivable assumptions for the sum will lead to meaningful or nontrivial solutions. 
Taking for instance R(v) = v 4 (A + v) or R(v) = v s (A + v)(B + v) with some unknown 
constants A and B leads in both cases to K2 = or j3 = 0, which we exclude for the above 
mentioned reasons. 

When demanding vanishing asymptotic boundary conditions for the u-field and its 



derivatives, the relations (4.9) and (4.11) imply that we satisfy the additional constraints 



2 2 

C^Kio Kirf ^2 /j \ 

lim u,ur,ucc = =>• Ki = — k 3 = — , hm v = . (4-16) 

C^ioo s ss 2c 2 2c c^±oo c 
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4.1.1 Type I solutions 

The simplest possible factorization for the sum in ( [4.13| ) is R(i 
up to the constraints 



«3 



(2f3a - 9a 2 - f3 2 ) 



K-2 



3V3c 2 {a-(3f 



16a 2 /3 ' '" z 16a 3 / 2 (-/3) 3 /2 ' 

c 3 (3a - P) (9a 2 - 6/3q + 5/3 2 ) 



and 



u(« — A) 4 , which holds 



v3c(/3 — a) 



(4.17) 



4a 3 / 2 ^=^ 



For a given specific model, i.e. fixed a, j3, 7, cp, this means that all remaining free parameters 
are fixed. Vanishing asymptotic boundary conditions for u require by the first two equations 
m 

CH) only one further relation f3 = —3a, despite the fact that one has to solve two 
constraining equations. Solving (|4.15| ) for the given factorization yields 



C - Co = ± 



aVx 



arc! anil ( 



A 



+ 



(4.18) 



First we focus on the real solution, which is obtained from (]4. 18f) for A > and < v < 
A. A further solution is obtained when replacing in ( 4.18| ) the arctanh by arcoth, which 
produces a real solution for < A < v. Again we are faced with the problem that we are 
not able to extract from (4.18) the function v(C)- Nonetheless, we can re-write ( [4.18 ) as 



v(C) = A + 



yMOVA 



arctanh 



MO 



±AV 2 ^/X(C-Co) 



(4.19) 



similarly as done in [43]. This is of course still not u(C), but this form is very useful 
to extract various types of information in an analytical manner. Simply by consider- 
ing the functions on the right hand side of (4.19|), we observe the asymptotic behaviour 



lim^ioo v(Q = A. In principle this is already sufficient to extract the qualitative be- 
haviour of i>(C)) since we also have assumed a simple factorization for the derivative of v 
in (4.13). We can, however, be more precise because the formulation ( 4.19| ) is also ideally 
suited to be solved numerically. In principle this possibility will become more powerful 
when considering more complicated scenarios. However, in these cases one encounters also 
more occurrences of the field v and it is not always obvious for which one the equation 
should be solved. 

For the purpose of a numerical study we discretise the equations as v n +i = F(v n ) 
and subsequently solve them iteratively for given values of (. The fixed points f/(C) of 
such discretised equations are known to be stable if and only if \F' [f/(C)]| < 1- We use 
this criterium to facilitate the numerical investigations. Concretely we solve the recursive 
equation 



v n+1 (() = A + 



arctanh 



i^/^C-Co) 



(4.20) 
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which converges very rapidly to a precision of ~ 10 5 typically already for less than 150 



iterations. We proceed similarly for the solution when < A < v. In figure 30 we depict 
two types of cusp solutions obtained in this manner. 
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Figure 30: Cusp solutions with asymptotically vanishing boundary conditions for the M-field of 
the Ito type equation with A = 3, c = 1, a = —1/3, j3 = 1, 7 = — 1/27 2 and Co = 0: (a) u-field; (b) 
M-field. 

In principle we could also proceed in this manner when taking complex initial condi- 
tions, but it is simpler to produce the contour plot in the way outlined at the beginning of 



section 3.1. For the same values of the parameters as in figure |30| we depict our results in 
figure [3l], indicating as before the imaginary parts of Co on some particular trajectories. 

As expected from ( |4.19| ) we observe in figures 31 that the complex solutions tend to 
the same asymptotic value as the real ones. The point A is an unstable focus in the 
"PT-symmetric and its broken version. In panel (b) more Riemann sheets are taken into 
account, revealing more substructure compared to panel (a). In the "PT-symmetric case we 
also observe the crucial feature that v*(Q = v(—Q and u*(Q = u (~ C)> which guarantees 
the reality of the energy as defined by the expression in (p..2|). 

4.1.2 Type II solutions 

Next we introduce an additional parameter B and assume the factorization of the form 
R{v) = v(v — A) 2 (v — B) 2 , which imposes the four constraining equations 



Ki 



K 2 



{p - 2a)c 2 {A 2 + AAB + B 2 ) a 
2o^ + 2 : 

2AB(A + B)a 3 
c(f3 - a) ' 



(4.21) 
(4.22) 



«3 



A + B 



2AB(A + B) 2 a 3 + c 3 (3a 



2/3) c (A 2 + \AB + B 2 ) 



c{a - 
\/3c(a 



6a 3 



2a 3 / 2 



(4.23) 
(4.24) 



This means our five constants K\, K2, K3, A, B are constrained by four equations, such that 
one of them remains free and thus allows us to adjust for some desired boundary conditions. 



- 35 - 



VT in complex nonlinear wave equations 




-2 024 -2 024 0246 

Rc v Rc i> Rc v 



Figure 31: Complex type I solutions with asymptotically vanishing boundary conditions for the 
w-field and corresponding w-held of the Ito type equation: (a) 'P7~-symmetric case with A = 3, 
c = 1, a = —1/3, (3 — 1 and 7 = — 1/27 2 ; (b) same values as in panel (a) for a single trajectory 
with Co = 0.04; (c) broken 'PT-symmetric case with A s» 3.054 — 0.783«, c = 1, a — —1/3, /3 = 1 — i 
and 7 = -1/27 2 - i/2. 



We exclude the trivial solutions A = B = A + B = 0as they all lead to K2 = 0. Solving 



(|4.15|) for this factorization gives 



C-Co 



(A-B)y/X 



arctanli ( -^7= 



B 



arctanh ( 
V VA 



(4.25) 



which reduces to ( f4.18| ) in the limit B — > A. Arguing as in the previous case, the real 
solutions are obtained from ( 4.1 8|) for A > and < v < A < B. For other configurations 
of the ordering we replace in ( 4,25| ) the arctanh by arcoth when the argument becomes 
greater than one. For the same reasons as in the previous section we isolate v from this 
equation and obtain 



v(C) = yltanh 2 



-VxVA(B -A)((- Co) + \l ^ arct anh 
2 \ B 



B 



(4.26) 



Simply by considering the functions on the right hand side of ( 4.26 ), we observe the asymp- 
totic behaviour lim^-i-oo w(C) = A. In this case the vanishing asymptotic boundary condi- 
tions for u can be implemented without additional constraints on the model defining pa- 
rameters a,/3,7,(/>. For other orderings of v,A,B we may also obtain lim^-i-oo v(C) = B. 
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When lim^ ±oov{() = A all constraints in ( |4.16|) and (|4.21| )- (|4.24D are satisfied with the 

choice 

c 2 (/3-3a) _ J n cyf^ 



A = -^——^= and B = --^_ (4.27) 



whereas when lim^_ s>= |- 00 v(Q = B we need to exchange A and B in ( 4.27| ). For the solution 



reported in [43| the possible values for a, /3,7 were restricted because the value of B was 
pre-selected, but as we have shown here this is not necessary. 



Choosing vanishing boundary conditions we depict the real solutions in figure 32 for 
the different regimes and real initial values Co- 
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Figure 32: Soliton, kink and cusp type solutions for the Ito type equations with a — —f3 = 
l/(2y/S), 7 = 1, c = —1, A = 2 and B = 4. The initial values are takes to be Co = — 10 solid (blue), 
Co = dashed (red) and Co = — 10 dotted (black), (a) -u-field; (b) u-field. 



Our findings disagree slightly with those in [43], where the case < v < \A\ < \B\ was 
reported to be of cusp type, whereas we observe that it is of a standard soliton nature. 
Our numerical findings are in agreement with the factorization of the right hand side of 
equation (4.13), which implies that (0) = and not infinity as needed for a cusp solution. 
We also find a new kink type solution in the region < \A\ < v < \B\ not reported by 
Kawamoto. In the regions < v < A < B and < A < B < v we observe explicitly the 
PT-symmetry £ — > —( + 2C , v — >• v and u — > u. 

Next we investigate some complex solutions by taking the initial values Co to be purely 
imaginary. For this type of the factorization the linearisation around the point A and B 
is straightforward as the square root in ( 4.15 ) can be taken. Parameterizing A = r^e A , 
B = r B e idB and A 
are computed to 



r\e x the eigenvalues of the Jacobian when linearized about v = A 



j k = ±y/r A r\ 

+i(-l) fc v /^ 



, 39 A e x \ e A e x 



(4.28) 



for k = 1,2. For the linearisation around v = B we obtain ( 4.28 ) with A <-> B. 

By tuning our free parameters we can produce any desired characteristic behaviour for 
the fixed points at A and B. For instance in the T^T-symmetric setting for A,B € R + we 
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s o 




Figure 33: Complex asymptotically constant type II trajectories for the "PT-symmetric Ito type 
system for purely complex initial values Co with a — —ft = 1/(2 v3), 7 = 1, c=— 1,^4 = 2 and 
B = 4. (a) u-field; (b) u-field. 



always obtain two real degenerate eigenvalues when A G 
v = A and v 



and therefore star nodes at 
B. In figure 33 we report an example of this type with j\ = j% = =f1/\/6- 



We also observe in figure 33 a that the complex trajectories surround the real solution 



1 



with the asymptotic point or points in common. For instance, the trajectory with Im £ c 
corresponds to a complexified version of a real soliton solution in the regime < v < A < B 
with asymptotic behaviour lim^-i-oo v (C) = A. In the n-plane the real solution becomes 
a cusp solution running off to infinity, whereas the complex solutions close. We may also 
identify complexified versions of the kink solutions in the regime < A < v < B, such 
as for instance the trajectory with ImCg = 5 in the v-plane with asymptotic behaviour 



lim, 



v(0 



A and lim, 



C- 



i u(C) = B. In the u-plane the qualitative behaviour 



remains the same with asymptotic behaviour lim ( ^^_ 00 u(C) = and lim^oo u(C) = —c/a. 
These features may also be derived analytically from ( }4.26| ). The reality of the energy is 
once more guaranteed be the symmetry u*(C) = V (~C) an d = u(— £). 

Similarly as for the trigonometric solution of the KdV equation, by changing the pa- 
rameters we can predict some periodic solutions. In the PT-symmetric setting this is 
achieved for A,B £ M + when taking A € K - , as in this case the two eigenvalues become 
purely complex and therefore the fixed point becomes a centre. For this type of solution 
this can be achieved either for the point v = A or v = B. We depict an example of such 
solutions in figure [34| with eigenvalues j'1/2 = ±il/\/6- In the v-plane we observe that 
the complex trajectories encircle the points A or B whereas in the u-plane the trajectories 
surround the points —c/a and for Im(" > or ImC < 0, respectively. 

Again we observe the symmetry relations f*(C) = v(— C) and = u(— C) ensuring 

the reality of the energy, but as for the trigonometric solution of the KdV equation we may 
compute the energy in this case explicitly. Taking the trajectories surrounding A in the 
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A o 




113 4 5 6 



1 

Pu=U 



Figure 34: Complex periodic trajectories for type II 'PT-symmetric Ito type system for purely 
complex initial values C with a — —(3 = l/(2\/3), 7 = — 1, c = — 1, A = 2 and B = 4. (a) u-field; 
(b) u-field. 





Figure 35: Complex periodic trajectories for type II broken 'PT-symmetric Ito type system for 
purely complex initial values Co with a = 1/(2^/3), ft = ^ 2l ^^~ 1 ^ 1 , 7 = i (— j + ^/— 1 — 1) 

c = -1, A = (1 + i)V2 and B = {2 + 2i)V2. (a) v-field; (b) u-field. 



v-plane to be the contour T we compute 

* I 



r /r V\yjv{v - A)(v - B) 



dv 



-IV- 



cA 2 + K 2 A + £ - + — 



ay/A(A - B 

For the trajectories surrounding i? in the u -plane we obtain in a similar way £r s 



(4.29) 
(4.30) 
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E Ta {A^B). 

Considering the amount of free parameters we have at our disposal, the expression for 



( 4.28 ) also suggests that we will be able to generate any type of fixed point even for the 
broken PT-symmetric scenario. Most unexpected is probably that we may even generate 
periodic orbits in that case. For this to happen we require 36 a + 6\ = tt and 6 a + 26 ] b+ 
6\ = it to hold. A solution to these equations, leading to the eigenvalues j\/2 = ±i2y/2, is 



presented in figure 35. 



As already seen for the complex KdV-system in section 2.1.2, we can also break the 
"PT-symmetry in a more controlled way and restore the reality of the energy given by 



( 4.30 ). Making for instance the parameter choice 



a = r n e 



A 



c(/3 - a) 
2 v / 3crV 2 ^/^ 



/3 = 3r a e l 4 : 
B = 2A, 



7 = (Ji- 



P—a 



a 



-1, (4.31) 
(4.32) 



with fjj £ R + and r a S M being unconstrained constants, we obtain for the energy of a 
periodic trajectory around the point A the real expression 



IT 

3rl 




5^ 
3" 



(4.33) 



Notice that in this case we have two free parameters available allowing to tune the real 
energy together with the model, despite the fact that any of these the Hamiltonians is 
neither Hermitian nor P7"-symmetric. We depict an example for such type of trajectory in 
figure |36l. We notice that the trajectories are qualitatively the same as those with complex 



energies depicted for instance in figure 35 




6 0.0 




0.4 
Re v 



0.8 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

Reu 



Figure 36: Complex periodic trajectories for type II broken PT-symmetric Ito type system with 
real energy Et a ~ —0.4275 for purely complex initial values C for the parameter choice in (4.31 



and (4.32) with r a = 2 and ^ = 2; (a) u-field; (b) it-field. 
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4.1.3 Type III solutions 

Next we assume that R(v) = v(v — A) 2 (B + v){C + v), which is achievable upon the 
constraints 



«3 



Q9 - 2a)c 2 + [A 2 + 2A(B + C) + £C]a 3 

yl[2£C + A(g + C)]a 3 
c((3 — a) 



(4.34) 



[c 2 - [A 2 + 2(B + C)A + BC]a 2 ]c 4(2A + B + C) [2BC + A(B + C)]a 3 f3c 3 



2a 2 



2(13 - a)c 



3a 3 ' 



A 



3V-^« 3 (a- P) 2 c 2 BC{B + C) 2 2BC 



Pa 3 (B + Cf 



B + C 



Now we have two free parameters at our disposal. In this case the solution of Q4.15D leads 
to 



c-c = ± 



1 



AyfBy/Xv 



F 



arcsm 



C 



+ n 



A 
B : 



arcsm 



C 



(4.35) 



with F [(j)\m] denoting an elliptic function of the first and II [n; (b\m\ an incomplete elliptic 
function. 

Choosing B and C conveniently we compute some trajectories similarly as in the 



previous subsection and depict our results in figure 37 



E 




Figure 37: Complex periodic trajectories for the type III 'PT-symnietric Ito type system for purely 
complex initial values £ with a — —(3 = 2^/6/19, 7 = — 1, c = —1, A = 1, B = 4 and C = 8. 



We identify some trajectories surrounding the point A in the u-plane for which we 
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compute the energy as 
E Ta 



7T- 



«[»«)]- = 



H[v] 



VX^(v - B)(v -C)v(v - A) 



dv 



y/(A-B)(A-C)A 



cA 2 + k 2 A + 



/3 



a a A/ 



(4.36) 
(4.37) 



As seen for the type II solutions we may also in this case break the 'P7~-symmetry and still 
obtain periodic solutions. Moreover, we can break the symmetry further, that is also on 
the level of the Hamiltonian, and still render the expression for Ej~ A real. Since in this case 
we have even an additional parameter C at our disposal it is conceivable that we might 
even have a free variable left in the expression for the real energy. We leave this question 
for a future investigation. 



An example for broken "P7~-symmetry is depicted in figure 38. 



A 




A 




Figure 38: Complex periodic trajectories for the type III with broken 'FT-symmetry for the Ito 
type system for purely complex initial values Co with a = —(3 = 2\/6/i9, 7 = —1 + i, c = — 1, 
A = 1, B = 4 and C = 8. 
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5. PT-symmetric deformations of the Ito type equations 

Since the Hamiltonian TLj admits the four different types of antilinear symmetries (|4.5|) - 
( |4,8| ), we have now the options to apply the deformation maps 5f and 57 introduced in 
( |3,1| ) to the fields u and v in different combinations. Accordingly we define four different 
PT-symmetric models with two deformation parameters e, [i suitable normalized by the 
following Hamiltonian densities 

n ++ = a _ d u z _ JV-U^y+i _ ^ {iVx y+i (5.x) 

e,n 2 g 1 + £ \ xj 1 + // x> V 1 



-Hi' = j^u(iv)^ - § u 3 - -^(iu.) 6 * 1 + %l (5.2) 



«5T = T^»(^)" +I - <m)2+e + i" 1 + 5* <5 - 4) 

By construction we have VTij ■ Ti £ 3 ^ ^ tGi,n with i,j G {+, — } and lim^-*.! Hf,^ = %i- 
The corresponding equations of motion resulting from ( |4.4| ) are 

%++ : + avv x + /Smis + ^u xxx , e = 0, TL+~ : u t + av fl v x + /Jun^ + ^u xxx>£ = 0, 
u t + a{uv) x + (fyv XXX)lx = 0, u t + "(w^ + ^xxx = 0, 



(5.5) 



W e ,» '■ u t + avv x + Pu e u x + ju xxx = 0, H E A( : u t + a^^ + /3n e n x + ju xxx = 0, 
w t + a(uv) x + (fyv XXX)lx = 0, u t + a(nu /J ) a; + = 0. 



(5.6) 



Naturally there exist also possibilities to construct non-Hamiltonian deformations. Noting 

for instance that the first equation related to is also invariant under VT '-\ we may 

combine it with the second equation resulting from W^Z and define the VT-\ -invariant 

system 

u t + avv x + f3uu x + ju xxx , £ = 0, ^ 
v t + a(uVfj) x + (pv xxx = 0. 
We have now numerous new physically feasible theories to be investigated. Here we will 
only present few examples. Having two parameters at our disposal allows to obtain more 
analytic expressions for the solutions. Technically we have two problems to overcome. First 
of all we have to decouple the equations and subsequently carry out all the integrations. 
In order to achieve the first aim we focus mainly on the case <j) = as this will allow to 
express u in terms of v in a simple manner. This choice will eliminate the deformation 
term in TLf^ and we therefore concentrate here on TLfZ to study the interplay between 
the two parameters. 

5.1 The model TLf~(a,f3,^f,u,v) 

As a further simplification we set the constant K2 to zero, which results in the first integra- 
tion of the equation for x>t- We may then integrate again and obtain a decoupled equation 
solely involving the field u 

.£+ 1 



7£ 



C 2 P s a 1 — u /a\j=% _2_ 
2 6 2 1 + jti Vc/ 



(5.8) 
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with the two fields related as 



a 



(5.9) 



To be able to carry out the final integration in (|5.8| ) we would like the right hand side 
to acquire the form of a factorizable polynomial in it. This is possible for several specific 
choices and with the additional choice of e we may even construct analytic solutions. 

This case is entirely reducible to the deformation of the KdV 7, u) when noticing 

that 



3a J 



2c 



■i £ -\u 



a 



3/2 



(5.10) 



where we used the relation between u and v as specified in (|5.9|) for fi = 1/3. 



1/2,1/2 



The choice \i = 1/2 converts the right hand side of ( p.8[ ) into a fourth order polynomial. 
Let us next assume the factorization is of the form \{u — A) 2 (u — B) 2 , which is indeed 
possible if the following constraints hold 

, 2 



A 



-in/4 a 4 

6c 3 7 



c b fi 6 - Uc^a 4 ? 



«3 



c 5 (c 2 /3 2 - 12a 4 )' 



384a 



12 



(5.11) 



, cV^v^vW 7 ^ c 3 /3a 4 T ^3 y^a 8 /? 2 - 8c 4 a 12 

A= 4a^ ' 5 = 4^ ' (5 - 12) 

Notice that all the free parameters are fixed in this case, as a consequence of having already 



pre-selected K2 = 0. We analyse this model in figure 39 



s 









Figure 39: Complex periodic trajectories for the VT ^ -symmetric system 1/2 f° r purely 

complex initial values £ with a = 1, j3 = 3,7 = 1, c= 1, A = (3 — \/3)/4 and B = (3 — v / 3)/4. 



We observe that the PT^ -symmetry manifests itself now through f*(C) = ~ V (~C) 

and u*(C) = u(— £). Furthermore we recognise that the asymptotic limits of the it-field are 
^4 and B. 
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Of course there are plenty more models one may explore. 



5.2 The non-Hamiltonian deformation (5.7) 



Let us briefly comment on one of the possibilities to construct deformations of non- 
Hamiltonian systems fl5.7|) . We consider again the case 4> = and K2 = and integrate the 
first equation in (|5.7|) for e = 1 and fi ^ 3 



c 2 P ^ a.c._2_l — it 3^ii 

K 3 + Kin + -u - -n J + -(-j^" 1 ^ u 1 -** 

2 b 2 a 3 — /u 



(5.13) 



The v-field is related to the w-field via (|5.9| ). The boundary conditions have to be treated 
by distinguishing different case. For instance, when (3 — fi) / (1 — fi) > the vanishing 
asymptotic boundary conditions for u and its derivatives demand that n\ = K3 = 0. 

The simplest case to consider is \i = — 1, which is just the KdV-system with a re-defined 
speed of the wave c — > c + a 2 /c. 

A further simple example is to take [x = 1/3 for which we can find an explicit solution 
u((). The right hand side of (|5.13| ) can be factorized into X(u — A) 2 (u — B) 2 with 

c 3 a 4 /? 3 - 9c 4 a 2 /3 (a 3 /3 2 - 9c 4 a) 2 , ^ 

A= 4aV Kl= 27c« Tm? ' (5 ' 14) 

, c 3 a 2 /3 - v^x/^a 4 /? 2 - 6c 10 a 2 „ a 2 /3c 3 + v^V^a 4 /? 2 - 6c 10 a 2 /rir , 
A= 3^ ' B = 3^ ' (5 - 15) 



for which (5.13) can be integrated further and solved for u. We find 



The system is easily linearized and qualitatively we find a similar behaviour as for the 
Hamiltonian systems, which we will, however, not present here in more detail. Thus, even 
though it is less clear how the 'PT-symmetry can be utilized, we find that they are not 
fundamentally different from the Hamiltonian systems. 

6. Conclusions 

The main focus of this paper was to investigate the effects of PT-symmetry and its breaking 
in complex nonlinear wave equations. In general we find that unlike as claimed for orbits 
in quantum mechanical one particle models |1C, 11, 13], trajectories in the complex 



plane of the field of a PT-symmetric nonlinear system are not fundamentally distinct from 
those associated to systems with spontaneously or completely broken 'PT-symmetry. Just 
from the general type of trajectory one can not conclude which type of setting one is 
considering. Of course invoking the information on the symmetry one can identify on this 
basis the PT-symmetric case over the broken ones. However, based on this criterium the 
spontaneously broken cases are indistinguishable from the completely broken ones. This 
behaviour extends to the type of fixed points one may find. We observed that essentially 
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all types of fixed points, except saddle points, may occur, irrespective of the PT-symmetry 
properties of the model. When identifying and varying certain constants of the model as 
bifurcation parameter the fixed point may undergo a Hopf bifurcation. 

It appears that the entire phase space is filled out in the broken case, thus suggesting 
a chaotic behaviour. However, from the Poincare-Bendixson theorem 1 we know that for a 
closed bounded and connected region in two dimensions this is impossible to occur. 

The nature of the fixed points also indicates that the models are not Hamiltonian in 
the real and imaginary part of u, as only saddle points and centres would emerge in that 
case, whereas we found different types of fixed points and the absence of saddle points. 

With regard to the energies we confirmed that fully PT-symmetric systems have real 
energies, which could be calculated explicitly in many cases. When breaking this symmetry 
spontaneously by complexifying some free parameters in the solutions we obtained complex 
energies. As expected, the models related to the complex conjugates of these parameters 
have complex conjugate energies. More surprising are the findings obtained in various 
models that one may regain the reality of the energy by breaking the symmetry further. 
For the type III solutions of the Ito type systems it is conceivable the one might have free 
parameters in the expressions for the energy even when fixing the model. We conjecture 
that these models possess a different kind of antilinear symmetry yet to be identified. 

We found that complex PT-symmetric soliton solutions behave similarly as their real 
counterparts, albeit in the complex plane. The one-solitons travel in the complex plane 
while maintaining their overall shape and the two-soliton solution can be associated to 
two separate one-soliton solutions in the distant past and future. However, when the PT- 
symmetry is broken the nature of the one-soliton solution changes into a breather, which 
only regains its shape after it has traveled a certain time and distance. The two-soliton 
solutions for this case can also be separated into these breathers in the distant past and 
future. The energy of the two-soliton solution was found to be the sum of the energies of 
the constituent one-soliton solutions in all scenarios. It would be interesting to investigate 
these features also for N-soliton solutions for N>2 and different types of systems, such as 
the Ito type. 

Besides guaranteeing the reality of the energy, PT-symmetry was noted to be useful for 
various other reasons. We found that the symmetry allows for a natural "e-prescription" 
to facilitate the computation of the energies in the complex plane by means of the Cauchy 
theorem. Identifying the different types of PT-symmetries it also allows to formulate new 
physically feasible models with real energies. This procedure constitutes a natural and 
more general framework for some models which have already been known before, such as 
the generalized KdV-equations with the modified one as a special case. 

Since many of the models discussed here are integrable, it is worth mentioning that 

most of the analysis carried out here for the Hamiltonian may also be performed for other 

conserved quantities of the same order in the fields having the same type of symmetry 

property. For some of the charges this is in fact not the case. Although in many deformed 

1 Poincare-Bendixson Theorem: Let ip t be a flow for a two dimensional dynamical systems and let 
V be a closed, bounded and connected set T> £M 2 , such that <fi t (T>) C T> for all time. Furthermore T> does 
not contain any fixed point. Then there exists at least one limit cycle in V. 
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cases the question of integrability has not yet been answered decisively, the statement also 
holds for the lowest charges which are certain to exist. 

The main difference between the deformed models and their undeformed counterparts 
is in general the occurrence of more and more Riemann sheets with increasing integer value 
for e. We showed that in many cases non-integer values and even negative values for e give 
rise to interesting solvable models. The overall structure of the trajectories and the nature 
of the fixed points is not fundamentally different, except that they usual extend over several 
Riemann sheets. Clearly we only presented here a limited number of solutions and many 
cases still need to be explored, especially for the second type of deformation for which even 
the factorization of the P(u)-function remains an open issue. 

When comparing the Ito type system with the KdV system we note that the former 
is not simply an add on to the latter, but gives rise to more complex structures. We found 
even for the undeformed case some new solutions such as those of kink type hitherto not 
reported. We also pointed out some minor discrepancies when compared to the literature. 
With regard to the deformations, which are all entirely new proposals, the conclusions are 
similar as for the deformed KdV-case. The major difference is that the interplay between 
the two deformation parameters allows for even more possibilities. Also in this case many 
possibilities remain unexplored. 

As indicated at the end of section 2.1.3 and 3.1.3 we can obtain simple quantum me- 
chanical systems as special cases from our analysis. For instance, the archetype deformation 
of the harmonic oscillator Hamiltonian H = p 2 + x 2 {ix) e can be obtained in various ways 
via the identification for the traveling wave u — > x and £ — > t. As discussed, the deformed 
models Ti~ yield precisely the quantum mechanical model with potential V(x) = x 2 (ix) £ , 
of which some cases were studied in |^5[ . From the models presented here there are more 
possibilities to arrive at such potential, as for instance also the non-Hamiltonian models 
Q5.7| ) give rise to such type of potentials as may be seen from ( |5.13| ). Exploiting these ob- 
servations allows to obtain many solutions and properties of these systems easily from our 
analysis overlooked up to now. Obviously, one may also construct new interesting quantum 
mechanical models in this way which have not been investigated so far. 

In addition, this opens up immediately the more general question of studying properties 
of these systems as continuous functions of e, rather than selecting just certain specific 
values as in this manuscript. Building for instance on the analogy with the potential 
systems should certainly reveal fundamentally different kinds of behaviour in some regions, 
such that for instance T~L~ will probably have a qualitatively different behaviour for negative 
values of e. We leave these type of questions for future investigations. 

Clearly it would be very interesting to extend these type of analysis to other non-linear 
field equations such as Burgers, Bussinesque, KP, generalized shallow water equations, 
extended KdV equations with compacton solution, etc. 
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A. The ten similarity classes for J 

For convenience we recall here the ten different similarity classes characterizing the fixed 
points of a two dimensional linear system ( |2.9[) by the eigenvalues ji and j% of the Jacobian 
matrix J at the fixed point. 



ji eR 


jl > 32 > 


unstable node 




32 < 3i < 


stable node 




32 < < ji 


saddle point 


ji = 32, diagonal J 


Ji > o 


unstable star node 




ji <o 


stable star node 


ji = j2, nondiagonal J 


ji >o 


unstable improper node 




ji < o 


stable improper node 


ji eC 


Reji > 


unstable focus 




Reji = 


centre 




Reji < 


stable focus 



Table 1: Nature of a fixed point as classified by the eigenvalues of the Jacobian. 
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